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The centrality and system size dependence of multiplicity fluctuations of charged particles pro- 
duced in nuclear collisions at 158 A GeV was studied by the NA49 collaboration. Centrality selected 
Pb+Pb collisions, semi-central C+C and Si+Si collisions as well as inelastic p+p interactions were 
analyzed. The number of projectile participants determined on an event-by-event basis was used to 
characterize the collision centrality. The scaled variance of the multiplicity distribution obtained in 
the forward rapidity region (1.1 < y c . m . < 2.6) shows a significant increase towards peripheral colli- 
sions. The results are similar for negatively and positively charged particles and about 50% larger 
for all charged particles. String-hadronic models of nuclear reactions without the fusion process do 
not reproduce the rise of fluctuations from central towards peripheral collisions. The measured cen- 
trality dependence can be reproduced in superposition models with the assumption of contributions 
from target participants to particle production in the forward hemisphere or in string models with 
fusion. 
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I. INTRODUCTION 

Nucleus-nucleus collisions at relativistic energies have 
been intensely studied over the last two decades. The 
main goal of these efforts is to understand the properties 
of strongly interacting matter under extreme conditions 
of high energy and baryon densities for which the cre- 
ation of a quark-gluon plasma (QGP) is expected [H, 0. 
Various collision characteristics and their collision en- 
ergy dependence suggest that a transient state of de- 
confined matter is created at collision energies as low 
as 30 A GeV @. Fluctuations of physical observables 
in heavy ion collisions have been a topic of interest for 
some years as they may provide important signals regard- 
ing the formation of a QGP. With the large number of 
particles produced in heavy ion collisions at CERN SPS 
and BNL RHIC energies it has now become feasible to 
study fluctuations on an event-by-event basis [4|. In a 
thcrmodynamical picture of the strongly interacting sys- 
tem formed in the collision, the experimentally studied 
fluctuations of particle multiplicities H, 0, HI, , mean 
transverse momenta [0, [H, HH , and other global ob- 
servables, are related to fundamental properties of the 
system, such as specific hea t [Til [l5| , chemical potential, 
and matter compressibility [16j . These, in turn, may re- 
veal information on the properties of the equation of state 
near the QCD phase boundary [TBL H~7L IT8j . 

The main objective of this work is to study how the 
multiplicity fluctuations change with the varying num- 
ber of nucleons participating in the collision, with cen- 
trality and system size. First results on the centrality 
dependence in Pb+Pb collisions at 158 A GeV published 
by the WA98 collaboration Q indicated that the scaled 
variance of the multiplicity distribution near midrapid- 
ity increases towards peripheral collisions. Our analysis 
extends such measurements to C and Si nuclei and care- 
fully addresses the contribution from fluctuations of the 
number of participants which can dominate the fluctua- 
tions of extensive quantities such as the produced particle 
multiplicity. Furthermore an extensive comparison of the 
results to available models is performed. The data sam- 
ple consists of centrality tagged minimum bias Pb+Pb, 
semi-central C+C and Si+Si as well as minimum bias 
p+p collisions registered by the NA49 detector at the 
CERN SPS. The multiplicity fluctuations were studied 
in the forward rapidity region (1.1 < y c . m . < 2.6). The 
collision centrality was characterized by the number of 
projectile participants derived from the projectile specta- 
tor energy measured in a forward calorimeter. Selection 
of narrow intervals in this energy minimizes the variation 
of the number of projectile participants. Although the 
number of projectile and target participants are closely 
correlated, the number of target participants cannot be 
tightly constrained with this procedure. The implications 
of this problem will be discussed in detail. 

The paper is organized as follows. In Sec. |TT] the 



method of measuring multiplicity fluctuations is intro- 
duced and discussed. The NA49 set-up is presented in 
Sec. IIIII Experimental procedures, in particular event 
and particle selection, detector acceptance and central- 
ity determination are discussed in Sec. IIVI The results 
on centrality and system size dependence of multiplicity 
fluctuations are presented in Sec. [V] In Sec. I VII our re- 
sults are compared with those of other experiments. Pos- 
sible explanations of our results are discussed in Sec. lVIII 
The paper closes with a summary and conclusions. 

II. MULTIPLICITY FLUCTUATIONS 
A. Observables 

Let P{N) be the multiplicity distribution, then 

(n) = J2np(n) (l) 

is the mean value of the distribution. The variance of the 
multiplicity distribution is defined as 

Var(iV) = J2( N - (N)) 2 P(N) = (N 2 ) - (N) 2 . (2) 

Note, that for a Poisson distribution the variance equals 
the mean value, Var (N) — (N) . Mean value and variance 
of multiplicity distributions are the only observables used 
in this analysis. 

B. Participants and spectators 

In a description of nuclear collisions the concept of 
participant and spectator nucleons is very useful. The 
participants are nucleons which are removed from the 
Fermi spheres of target or projectile nuclei due to the 
collision. The remaining nucleons are called spectators. 
In the case of central nucleus-nucleus collisions, where the 
impact parameter b is relatively small, almost all nucle- 
ons participate in the collision. In particular, the number 
of projectile participants, Np ROJ , approximately equals 
the total number of projectile nucleons Np R ~ A. 
If the collision is peripheral (with large impact param- 
eter b) almost all nucleons are spectators, N PROJ <C A. 
The number of projectile spectators N RR ^ is given by 

N PROJ _ A _ atPROJ 

C. Multiplicity fluctuations in superposition 
models 

Superposition models of nuclear collisions, which are 
frequently used, assume that secondary particles are 
emitted by independent sources. A prominent example 
is the Wounded Nucleon Model [lj| , in which the sources 
are wounded nucleons, i.e. the nucleons that have inter- 
acted at least once. In this model the number of partici- 
pants would be equal to the number of wounded nucleons. 
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In superposition models the total multiplicity is given 



N s 



(3) 



where Ns denotes the number of sources and rm is the 
multiplicity of secondaries from the i— th source. When 
the sources are identical and independent from each 
other, the mean total multiplicity equals 



(N) = (N s )(m) 



(4) 



where (Ns) is the mean number of sources and (m) is 
the mean multiplicity from a single source. The variance 
of the multiplicity distribution is 



Var(iV) = (ATs)Var(m) + (m) 2 Var(^V s ), 



(5) 



where Var(m) and Vai(Ns) denote the variance of the 
distribution of single-source multiplicity and the variance 
of the distribution of source number, respectively. 

The scaled variance of the multiplicity distribution 
Vai(N)/ (N) is a useful measure of the multiplicity fluc- 
tuations. From Eqs. (HJ and ([5]) one gets 



Var(iV) Var(m) 
(N) = (m) 



Var(jVs) 
' (N S ) ■ 



(6) 



Thus, in superposition models the measured scaled vari- 
ance of the multiplicity distribution is the sum of two 
contributions. The first one describes the multiplicity 
fluctuations from a single source while the second one 
accounts for the fluctuations of the number of sources. 
To infer the fluctuations of main interest given by the 
first term in Eq. ^ we try to minimize the contribution 
of the second term. In this experiment events with a 
fixed number of projectile participants can be selected in 
which, however, the number of target participants can 
still fluctuate. The effect of the fluctuation of target 
participant number at fixed number of projectile par- 
ticipants was theoretically studied in [2(j. Multiplicity 
fluctuations are measured in a forward rapidity window 
in which one mainly expects particles produced from pro- 
jectile participants. Nevertheless the effect of remaining 
fluctuations of the number of target participants can re- 
main important, if participants produce particles over a 
wide rapidity range. We return to this point in Secs. lIVFl 
andEIIl 



III. NA49 EXPERIMENTAL SET-UP 

The NA49 experiment is a large acceptance hadron 
spectrometer at the CERN-SPS (see Fig. Q] and (U) 
used to study the final hadronic states produced in col- 
lisions of beam particles (Pb directly from the SPS, C 
and Si via fragmentation of the primary Pb beam, p as 
secondary beam from the 450 GeV proton beam in the 



SPS) with a variety of fixed targets. The main tracking 
devices are four large volume Time Projection Chambers 
(TPCs) which are capable of detecting 60% of some 1500 
charged particles created in a central Pb+Pb collision at 
158 A GeV. Two of them, the Vertex TPCs (VTPC-1 
and VTPC-2), are located inside the magnetic field of 
two super-conducting dipole magnets (1.5 and 1.1 T, re- 
spectively) and two others (MTPC-L and MTPC-R) are 
positioned downstream of the magnets symmetrically to 
the beam line. The results presented here are analysed 
with a global tracking scheme [12] , which combines track 
segments that belong to the same physical particle but 
were detected in different TPCs. The NA49 TPCs allow 
precise measurements of particle momenta p with a res- 
olution of a(p)/p 2 ~ (0.3 - 7) • 10~ 4 (GeV/c)" 1 . The 
set-up is supplemented by two Time of Flight (TOF) 
detector arrays and a set of calorimeters. 

The targets - C (561 mg/cm 2 ), Si (1170 mg/cm 2 ) discs 
and a Pb (224 mg/cm 2 ) foil for ion collisions and a liq- 
uid hydrogen cylinder (length 20.29 cm) for elementary 
interactions - are positioned about 80 cm upstream from 
the VTPC-1. 

Pb beam particles are identified by means of their 
charge as seen by a Helium Gas-Cherenkov counter (S2') 
in front of the target. The protons were identified by 
Cherenkov counters farther upstream of the target. The 
study of C+C and Si+Si reactions is possible through 
the generation of a secondary fragmentation beam which 
is produced by a primary target (1 cm carbon) in the 
extracted Pb-beam. Setting the beam line momentum 
accordingly, a large fraction of all Zj A — 1/2 fragments 
are transported to the NA49 experiment. On-line selec- 
tion based on a pulse height measurement in a scintillator 
beam counter (S2) is used to select particles with Z = 6, 7 
(C, N) and Z = 13, 14, 15 (Al, Si, P). Off-line clean-up is 
achieved by using in addition the energy loss measured 
by beam position detectors (BPD- 1/2/3 in Fig. [TJ) . These 
detectors consist of pairs of proportional chambers and 
are placed along the beam line for a precise measurement 
of the transverse positions of the incoming beam parti- 
cles. 

For p beams, interactions in the target are selected 
by anti-coincidence of the incoming beam particle with a 
small scintillation counter (S4) placed on the beam tra- 
jectory between the two vertex magnets. For p+p inter- 
actions at 158 A GeV this counter selects a (trigger) cross 
section of 28.23 mb out of 31.78 mb of the total inelastic 
cross section [23[ . For Pb-ion beams an interaction trig- 
ger is provided by anti-coincidence with a Helium Gas- 
Cherenkov counter (S3) directly behind the target. The 
S3 counter is used to select minimum bias collisions by 
requiring a reduction of the Cherenkov signal by a factor 
of about 6. Since the Cherenkov signal is proportional 
to Z 2 , this requirement ensures that the Pb projectile 
has interacted, with a minimal constraint on the type of 
interaction. This set-up limits the triggers on non-target 
interactions to rare beam-gas collisions, the fraction of 
which proved to be small after cuts, even in the case of 
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FIG. 1: The experimental set-up of the NA49 experiment with different beam definitions and target arrangements. 



TABLE I: The number of events and the fraction of the total 
inelastic cross section selected by the on-line trigger for data 



sets used in this analysis. 

Data Set No of events a/a inel 

p+p 319 000 0.9 

C+C 51 000 0.153 

Si+Si 59 000 0.122 

Pb+Pb 165 000 0.6 



peripheral Pb+Pb collisions. 

The centrality of A+A collisions is selected by using in- 
formation from a Veto Calorimeter (VCAL), which mea- 
sures the energy of the projectile spectator nucleons. The 
geometrical acceptance of the Veto Calorimeter is ad- 
justed in order to cover the projectile spectator region 
by a proper setting of the collimator (COLL). 

Details of the NA49 detector set-up and performance 
of tracking software are described in [2l[ . 
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FIG. 2: Rapidity distribution of negatively charged particles 
produced in central Pb+Pb collisions at 158 A GeV within the 
NA49 acceptance (full line) in comparison with an acceptance 
corrected distribution from [241 ] (dotted line). 



IV. DATA SELECTION AND ANALYSIS 

A. Data sets 

Multiplicity fluctuations are studied for negatively, 
positively and all charged particles selecting events 
within narrow intervals of energy measured by the VCAL 
(predominantly energy of projectile spectators). The 
experimental material used for the analysis consists of 
samples of p+p, C+C, Si+Si and Pb+Pb collisions at 
158 A GeV. The number of events in each sample is given 
in Tab. HI For Pb+Pb interactions a minimum bias trig- 
ger was used allowing a study of centrality dependence. 



B. NA49 acceptance 

The NA49 apparatus detects mainly the particles pro- 
duced in the forward rapidity hemisphere. Fig. [3 shows 
the rapidity distribution of negatively charged particles 
in the NA49 acceptance calculated assuming that all 
particles are pions. This distribution is compared to 
an acceptance corrected rapidity distribution of nega- 
tively charged particles from [2J|. As seen, the rapid- 
ity region where almost all particles are measured is 
1-1 < Vcm, < 2.6. 

The NA49 detector was designed for a large acceptance 
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in the forward hemisphere. However, also in this region 
the geometrical acceptance is not complete. The accep- 
tance limits in the transverse momentum and azimuthal 
angle were parameterized by a simple function 



Pt{4>) 



1 



A+% 



B, 



(J) 



where the values of A, B and C depend on the rapid- 
ity interval, see [Toj |. Only particles within the curves 
are used in this analysis. This well defined acceptance 
is essential for later comparison of the results with mod- 
els and other experiments. Only forward rapidity tracks 
(1.1 < Dc.m. < 2-6, rapidity calculated assuming pion 
mass for all particles) with 0.005 < pr < 1.5 GeV/c have 
been used in this analysis. 



C. Event and particle selection 

The aim of the event selection criteria is to reduce con- 
tamination from non-target collisions. The primary ver- 
tex was reconstructed by fitting the intersection point of 
the measured particle trajectories. Only events with a 
proper quality and position of the reconstructed vertex 
were accepted for further analysis. The vertex coordinate 
Z along the beam had to satisfy \Z— Zq\ < AZ, where the 
nominal vertex position Zq and cut parameter AZ val- 
ues are: —579.5 and 5.5 cm, —579.5 and 1.5 cm, —579.5 
and 0.8 cm, -581.2 and 0.6 cm for p+p, C+C, Si+Si 
and Pb+Pb collisions, respectively. The vertex position 
in the transverse X, Y coordinates had to agree with the 
incoming beam position as measured by the BPD detec- 
tors. 

In order to reduce the contamination of particles from 
secondary interactions, weak decays and other sources of 
non-vertex tracks, several track cuts were applied. The 
accepted particles were required to have measured points 
in at least one of the Vertex TPCs. A cut on the extrap- 
olated distance of closest approach to the fitted vertex of 
the particle at the vertex plane was applied (\dx\ < 4 cm 
and \dy\ < 2 cm). Moreover the particle was accepted 
only when the potential number of points (calculated on 
the basis of the geometry of the track) in the detector ex- 
ceeded 30. The ratio of the number of points on a track 
to the potential number of points had to be higher than 
0.5 to avoid split tracks (double counting). 



D. Centrality selection 

In order to reduce the effect of fluctuations in the num- 
ber of participants (and thus of particle sources N$), the 
multiplicity fluctuations were analyzed in narrow central- 
ity bins defined by the energy measured in the VCAL. 
This procedure minimizes the variation of the number of 
projectile participants Np ROJ . Although the number of 
projectile and target participants are closely correlated, 



the number of target participants cannot be constrained 
with the NA49 detector. This leads to a remaining fluc- 
tuation in the total number of participants, i.e. particle 
sources N§. 

For C+C and Si+Si interactions, three and five narrow 
centrality bins were selected respectively. In the case of 
Pb+Pb collisions, 58 narrow centrality bins were chosen. 
For each centrality bin the number of projectile partici- 
pants Np ROJ was estimated by 



N PROJ 



A 



E 



Veto 



Elab 



(8) 



where Eveto is the energy deposited in the Veto 
Calorimeter; Elab is the energy carried by single pro- 
jectile nucleons. The resolution of Eyeto, which is re- 
calculated into the resolution of N RROJ , is discussed in 
Sec. EH 

The effect of smearing of energy carried by spectator 
nucleons due to their Fermi momentum was estimated. 
The effect was found to be significantly smaller than the 
resolution of the VCAL and for this reason can be ne- 
glected. 

We also studied the contribution of non-spectator par- 
ticles to the energy measured by the veto calorimeter [43[ 
with detailed simulations using VENUS and UrQMD 
events. The net result was a reduction of the scaled 
variance of the multiplicity distribution by about 0.3 for 
central Pb+Pb collisions when the non-spectator contri- 
bution to the VCAL energy is taken into account. With 
increasing number of spectators the effect becomes much 
smaller due to this effect. No corrections were applied 
to the experimental results, nor were the non-spectator 
contributions taken into account in the calculations of 
model predictions. 



E. Multiplicity distributions 

The multiplicity distribution depends on the selected 
Eveto interval (its position Ey and width AEy ; see Fig. [3] 
for definitions) and the kinematic acceptance selected for 
the analysis. 

In the centrality intervals and acceptance selected for 
this analysis, multiplicity distributions show approxi- 
mately Poissonian behavior for p+p and central Pb+Pb 
collisions (see Fig. @|. For semi-peripheral collisions the 
multiplicity distribution is significantly broader than a 
Poissonian distribution. 

In Figs. [5] and [S] the measured mean value and the vari- 
ance of the multiplicity distributions as a function of the 
number of projectile participants Np ROJ are presented. 
The quantities are not corrected for the fluctuations of 
the number of projectile participants, as discussed in the 
next subsection. One can see from these plots that the 
mean multiplicity shows an approximately linear depen- 
dence on the number of projectile participants, whereas 
the variance of the multiplicity distributions calculated 
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FIG. 3: (color online) Distribution of energy deposited in 
the Veto Calorimeter for minimum bias Pb+Pb collisions at 
158 A GeV. An example of an Eyeto interval is shown; the 
interval is determined by its central value Ey and the width 
AEy. The width of the histogram bins corresponds to the 
width of the AEy intervals used for the fluctuation analysis. 
The upper horizontal scale shows the corresponding values of 
the number of projectile participants Np ROJ . 



from the data exceeds the variance of the Poisson distri- 
butions. 



F. Scaled variance of the multiplicity distribution 

The scaled variance of the multiplicity distribution de- 
pends, among other variables, on the width of the energy 
interval AEy selected in the VCAL. For very broad AEy 
intervals the measured scaled variance of the multiplicity 
distribution has a large value because of significant fluc- 
tuations in the number of projectile participants N PROJ . 
Narrowing the AEy interval results in decreasing fluctu- 
ations in the number of projectile participants and conse- 
quently in a reduction of the scaled variance as shown in 
Fig. [3 For AEy smaller than about 1 TeV the measured 
scaled variance of the multiplicity distribution is almost 
independent of AEy. Note that even for very small val- 
ues of AEy the number of projectile spectators fluctuates 
due to the finite resolution of the VCAL. In the follow- 
ing the scaled variance will be calculated in the Ey eto 
intervals AEy — 500 GeV (the choice of the width of the 
interval is a compromise between minimizing correction 
for the finite interval width and sufficient statistics) and 
then corrected for fluctuations in the number of projectile 



participants due to the finite width of the Ey eto interval 
and the finite resolution of the Veto Calorimeter. 

Within the superposition model leading to Eq. ([6]), the 
correction 5 to the scaled variance, which takes into ac- 
count the finite value of AEy and the calorimeter reso- 
lution, is calculated as 



(AO (Var A (E Veto ) + Var fl (E Veto )) 



E 



BEAM 



— (Eyeto)) 



(9) 



where Var a {Ey eto ) is the variance of Ey eto due to the 
finite width of the Ey eto bin, Var^ (Ey eto ) is the vari- 
ance of Ey eto due to the VCAL resolution, (Ey eto ) is the 
mean value of Ey eto in the bin and Ebeam = 158 A GeV 
is the total beam energy. Var a (Ey e to) was calculated 
from the distribution of Ey eto energy in a given interval. 
V&r R (E Ve to) = cr 2 (E Veto ) with E Veto expressed in GeV 
was parameterized using calibration measurements [25[ 
as: 



C {Eyeto) 
Eyeto 



2.85 



1G 



VE- 



Veto 



Ey e 



(10) 



The resolution of Ey e to is easily propagated to the reso- 
lution oiN^ ROJ using Eq. ©. 

Finally, the corrected scaled variance is calculated as 



Var (n) _ Var (AT) 



(11) 



where Var (N) is the measured variance, (N) the mea- 
sured mean value of the multiplicity distribution in a 
given Eyeto bin, and 6 represents the corrections for fluc- 
tuations in the number of projectile participants. The 
dependence of 5 on AEy for some values of Ey are pre- 
sented in Fig. [5J while the total correction S as a function 
of number of projectile participants Np ROJ is presented 
in Fig. [9] Using the same procedure, corresponding cor- 
rections were determined for C+C and Si+Si collisions. 
The corrected values of the scaled variance from Eq. (fTTj) 
are plotted in Fig. [TTJ 

Let us discuss in more detail the physical meaning 
of the scaled variance displayed in Fig. QT] which is the 
main result of our study. As already stressed, fluctua- 
tions of the number of projectile participants contribute 
to the multiplicity fluctuations of interest. However, the 
scaled variance of the multiplicity distribution can be 
split within the superposition model according to for- 
mula ([5]), that is into the contributions coming from a 
single participant (or equivalently from a fixed number 
of participants) and from the varying number of partici- 
pants - the first and the second term in Eq. ([6]), respec- 
tively. Since the correction S represents the contribution 
from the fluctuating number of projectile participants, 
the corrected scaled variance shown in Fig. [11] describes 
fluctuations at fixed number of projectile participants. 

If the number of projectile participants uniquely deter- 
mines the number of particle sources, which contribute 
to the rapidity interval under study, the scaled variance 
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FIG. 4: (color online) Multiplicity distributions of negatively charged particles obtained for AEv = 500 GeV for central Pb+Pb 
collisions with number of projectile participants Np HOJ = 178 (left panel); semi-peripheral Pb+Pb collisions with Np ROJ — 39 
(mid-panel) and p+p interactions with Np HOJ — 1 (right panel). 
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FIG. 5: (color online) The uncorrected mean value of the 
multiplicity distribution for negatively charged particles as a 
function of the number of projectile participants, (rapidity 
interval 1.1 < y c . m . < 2.6) 



FIG. 6: (color online) The uncorrected variance of the multi- 
plicity distribution for negatively charged particles as a func- 
tion of the number of projectile participants. 



from Fig. [TTJ represents the multiplicity fluctuations of 
particles coming from a single source. However, the num- 
ber of particle sources can fluctuate even at fixed num- 
ber of projectile participants. For example, the number 
of strings, which decay into particles, is a random vari- 
able in some models. The number of particle sources can 
also fluctuate due to the fluctuating number of target 
participants which are not observed in the NA49 detec- 
tor. Thus, we stress that the scaled variance shown in 
Fig. [TTJ represents the multiplicity fluctuations at fixed 
number of projectile participants and may still contain 



some contribution from fluctuations in the number of 
particle sources. 



G. Statistical errors 

In general, the variance of a quantity expressed by the 
fraction u> = x/y is given as 



Var(w) 



(y) 



2 r 



y 

(y) 2 



2Cov(x, y) 



(12) 
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FIG. 7: Measured scaled variance of multiplicity distributions 
of negatively charged particles for minimum bias Pb+Pb col- 
lisions as a function of the interval width AEv of the selected 
VCAL energy for various positions Ev of the interval. The 
vertical line shows the width of the interval, which was used 
for further analysis. See text for details. 
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FIG. 8: The corrections for Veto Calorimeter resolution and 
finite interval width as a function of interval width AEv for 
various positions Ev of the interval. 



where the symbol (. . .) means an average over events, a 
denotes the statistical error of x or y, and Cov(x, y) is 
the covariance between x and y. The statistical error of 
the quantity lu is given by a (oj) — yVar (x/y). 

For the calculation of the statistical errors of the scaled 
variance, we have x = Var (N), y = (N) and 



Var (AH r- 

a* = Jj- V 2 + 72, 

Cov (x, y) 



M3 (N) 



(13) 



where 72 = /U4 (N) /fx^ (N) — 3 is the kurtosis, fik (N) is 
the k— th central moment of the multiplicity distribution, 
and L is the number of events. 

The statistical errors of our results are usually much 
smaller than the systematic uncertainties discussed in 
Sec.UvTl 



H. Stability for event cuts and track selection 

In this section the stability of our results with respect 
to variations of the event and track selection criteria is de- 
scribed. All tests were performed for negatively charged 
particles produced in Pb+Pb minimum bias collisions for 
the centrality bin corresponding to the number of projec- 
tile participants N^ ROJ = 39. 

The event and track selection criteria are designed to 
reduce the contamination from the background. To check 
the stability of the results obtained for Var (n neg ) / {n neg ) 




FIG. 9: The total correction applied to the measured scaled 
variance of negatively charged particles as a function of the 
number of projectile participants Np ROJ . 



the cut parameters were varied within reasonable ranges. 
In addition, results obtained from the analysis of data 
taken at two different magnetic field polarities as well as 
during different running periods have been compared. 

Event cuts are used to reject contamination by non- 
target interactions. Fig. [TUh shows the scaled variance 
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FIG. 10: Scaled variance of the multiplicity distribution of 
negatively charged particles for the centrality bin correspond- 
ing to the number of projectile participants Np ROJ = 39 as 
a function of: (a) maximally allowed difference AZ between 
fitted main vertex and target position; (b) maximally allowed 
difference between the transverse beam position measured by 
the BPD and the fitted main vertex position (\Xbpd — Xfit\) 
(simultaneously the difference in \Ybpd — Yfit\ was required 
to be below 2/3 of the upper limit in \Xbpd — Xfit\); (c) 
maximally allowed distance between the reconstructed pri- 
mary vertex and the track in the target plane \dx I (simultane- 
ously the deviation in \dy \ was required to be below 0.5- \dx I); 
(d) minimum value of the ratio of the number of points on 
the track Np to the potential number of points Nmp- The 
vertical lines indicate the cuts used to obtain the results. 



as a function of the maximally accepted distance AZ 
between fitted and nominal Z position of the primary 
vertex. As one can see, the scaled variance of the multi- 
plicity distribution is stable with respect to a small con- 
tamination by non- vertex tracks. 

Fig. [TUb shows the scaled variance as a function of the 
maximally allowed difference of the X and Y position of 
the fitted main vertex from the transverse beam position 
given by the beam position detectors (cf. Fig.[T]). Also in 
this case the scaled variance is stable. 

The majority of tracks selected by the track selection 
criteria are main vertex tracks and the remaining frac- 
tion (ps 10%) originates predominantly from weak decays 
and secondary interactions with the material of the de- 
tector. To estimate the influence of this contamination 
on the multiplicity fluctuations, the maximally accepted 
distance between the extrapolated track and the recon- 
structed primary vertex in the target plane was varied 
(see Fig.UDfc). 



Losses of tracks due to reconstruction inefficiency and 
track selection cuts influence the measured multiplicity 
fluctuations. In order to estimate this effect, the depen- 
dence of Var (n neg ) / (n neg ) on the ratio of the number of 
points on a track to the potential number of points was 
determined (see Fig. ITOtO. 

In summary, the values of the scaled variance appear 
stable with respect to reasonable variations in the event 
and track selection cuts. 



I. Systematic uncertainties 

There are several sources of systematic uncertainties of 
our results. The systematic error due to the contamina- 
tion of non-vertex interactions, tracks from weak decays 
and secondary interactions as well as reconstruction inef- 
ficiencies and biases were estimated by varying event and 
track selection cuts and simulations as discussed in the 
previous section. The systematic error due to the above 
mentioned effects was estimated to be 10%. 

The main source of a possible systematic bias, how- 
ever, is the uncertainty in the measurement of the energy 
deposited in VCAL. To estimate the influence of this un- 
certainty on the scaled variance of the multiplicity dis- 
tribution the analysis of so-called beam only events was 
performed. Here, the Pb beam hits into VCAL with- 
out interactions with target nuclei. The measured res- 
olution was a (Eveto) /Eveto = 2.5/ -JE Veto (E Veto in 
GeV) [25j. However, the beam impinges on a small spot 
in the VCAL, and the resolution estimated by using beam 
only events does not take into account the VCAL non- 
uniformity To estimate it, the correlation between 
signals deposited in different sectors of the VCAL were 
studied. From the result the resolution (see Eq. (flT)|) ) 
of the Veto calorimeter was derived [25| . The system- 
atic error due to the finite resolution and non-uniformity 
of the VCAL was estimated to be smaller than 10% for 
N PRO.J > 20 and 25% f or N PROJ < 20. 

Finally, the systematic error of the results in Pb+Pb 
collisions was set to be equal to 10% for N PROJ > 20 and 
25% for N PROJ < 20. The resulting estimates are shown 
by the outer error bars in Fig. 1111 For semi-central C+C 
and Si+Si collisions we assume 10% systematic error (cf. 
Figs. [L21TT5)) . In p+p interactions the centrality selection 
does not apply and the systematic error was estimated 
by varying event and track selection cuts yielding 10% 
error. 



V. RESULTS 

The results discussed in this section refer to the ac- 
cepted particles, i.e. particles that are registered by 
the detector and pass all kinematic cuts and track se- 
lection criteria. The data cover a broad range in px 
(0.005 < px < 1-5 GeV/c). The center-of-mass ra- 
pidity of accepted particles is restricted to the interval 



10 



1.1 < y c . m . < 2.6 where the azimuthal acceptance given 
by Eq. ([7]) is large. 



A. Multiplicity fluctuations in Pb+Pb collisions 

The corrected scaled variance of the multiplicity distri- 
bution for negatively, positively and all charged accepted 
particles produced in minimum bias Pb+Pb collisions as 
a function of centrality is shown in F ig. [TT1 (cf. [!, [261') 
and compared to HIJING [27[, HSD [il]7~UrQMD (29| 
and VENUS [3(| simulations. In the model calculations 
a realistic simulation of the determination of the number 
of projectile participants was included. The models pro- 
duce approximately Poissonian multiplicity distributions 
independent of centrality. The data points, in contrast, 
indicate a strong increase towards peripheral collisions 
possibly with a maximum at about Np ROJ ~ 30 projec- 
tile participants. We note that by chance the scaled vari- 
ance is approximately 1 for p+p collisions at 158 A GeV, 
although the multiplicity distribution in p+p collisions 
is not Poissonian both at lower and higher energies (cf. 
Sec. EJ and [HI). 

The scaled variances for positively and negatively 
charged particles are similar. The corresponding values 
for all charged particles are larger. Assuming that nega- 
tively and positively charged particles in the experimen- 
tal acceptance are correlated with a correlation factor p, 
one gets 



Ya,r(n ch ) Var(n neg ) 



{rich) 



(1+P) 



Var (n 



pos J 



L neg / 



L pos I 



(14) 



where n neg , n pos and n c h are multiplicities of negatively, 



positively and all charged particles, respectively, and, for 
simplicity, we assumed (n neg ) = (n pos ). As seen, the 
scaled variance for charged particles is 1 + p times as 
large as that for positives or negatives (see also [32]|). 
The positive value of the correlation factor required by 
the data (p « 0.5) may reflect the effect of electric charge 
conservation. 
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FIG. 11: (color online) The scaled variance of the multiplicity 
distribution for negatively (upper panel), positively (middle 
panel) and all (bottom panel) charged particles as a function 
of the number of projectile participants Np ROJ compared 
with model simulations in the NA49 acceptance (HSD and 
UrQMD predictions were taken from [3^|). The statistical 
errors are smaller than the symbols (except for the most pe- 
ripheral points). The horizontal bars indicate the systematic 
uncertainties (the upper one for the second point is off scale). 



B. Multiplicity fluctuations in C+C and Si+Si 
collisions 

The corrected scaled variances of multiplicity distribu- 
tions for negatively, positively and all charged accepted 
particles produced in semi-central C+C and Si+Si colli- 
sions as a function of centrality are presented in Figs. 1121 
and 1131 The data are compared with results from a HI- 
JING simulation. The scaled variance for negative and 
positive hadrons shows almost Poissonian fluctuations in 
C+C as well as in Si+Si collisions, similarly to the ob- 
served multiplicity fluctuations in central Pb+Pb inter- 
actions. The multiplicity fluctuations for all charged ac- 
cepted particles are much larger than for the like sign, 
as also observed in Pb+Pb collisions. Since data on pe- 
ripheral C+C and Si+Si collisions are not available, it 



is difficult to speculate about the centrality dependence 
of the scaled variance. However, it seems that there is a 
scaling seen in the variable / = Np ROJ /A, the fraction of 
nucleons participating in the collision (Fig. 1 14)) . In order 
to validate this very suggestive scaling data on minimum 
bias collisions of C+C and Si+Si would be necessary. 



VI. COMPARISON WITH RESULTS FROM 
OTHER EXPERIMENTS 

Assuming that the produced particles are emitted in- 
dependently in momentum space, the mean multiplicity 
in a limited acceptance (N) can be expressed through 
the mean multiplicity in the full acceptance (N p —i) as 
(N) = p(Np—i), where p is the fraction of registered 
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FIG. 12: (color online) The scaled variance of the multiplic- 
ity distribution of negatively (upper panel) , positively (middle 
panel) and all (bottom panel) charged particles produced in 
semi-central C+C collisions as a function of the number of 
projectile participants N PROJ compared with HIJING simu- 
lations in the NA49 acceptance. The horizontal bars show the 
systematic uncertainties while the statistical errors are repre- 
sented by the vertical bars. Except for the most peripheral 
points the statistical errors are smaller than the symbols. 



FIG. 13: (color online) The scaled variance of the multiplic- 
ity distribution of negatively (upper panel) , positively (middle 
panel) and all (bottom panel) charged particles produced in 
semi-central Si+Si collisions as a function of the number of 
projectile participants N PROJ compared with HIJING simu- 
lation in the NA49 acceptance. The horizontal bars show the 
systematic uncertainties while the statistical errors are repre- 
sented by the vertical bars. Except for the most peripheral 
points the statistical errors are smaller than the symbols. 



particles. The variance of the multiplicity distribution 
in a limited acceptance can be written as Var (N) = 
p 2 Var(iVp = i) + (N p=1 )p(l -p). The probability that a 
given fraction of produced particles is registered is given 
by the binomial distribution. Using the above formulas, 
one easily finds the scaled variance in a limited accep- 
tance 



Var (N) 



1 



/ Var (JV p=1 ) 



i 



(15) 



which depends linearly on the fraction p of registered 
particles. 

In order to compare to results obtained by other ex- 
periments, the NA49 results were extrapolated to the 
full azimuthal acceptance in the rapidity interval 1.1 < 
Uc.m. < 2.6. The extrapolation was done by means of 



Eq. (jT5|) and values of p obtained from a HIJING simula- 
tion which included the NA49 acceptance filter. Fig. [15] 
shows the scaled variance of the multiplicity distribution 
of all charged particles for N PROJ = 29 from the NA49 
and WA98 @ experiments plotted versus the fraction p 
of registered particles in the NA49 experiment and in the 
WA98 experiment [j| . The latter measured the multiplic- 
ity fluctuations of charged particles produced in Pb+Pb 
collisions at 158^4 GeV in the central rapidity region. 
The predicted linear dependence of Eq. (fT5|) is seen to 
agree with WA98 results [j| . The WA98 experiment used 
a centrality determination method which is different from 
NA49, i.e. the total transverse energy Et in the forward 
rapidity region was used to fix the centrality in the WA98 
experiment. Nevertheless the results of NA49 and WA98 
at the same value of p agree (see Fig. IT51) . 



12 



ra 1.5 
> 

1 

0.5 



2 
1.5 
1 

0.5 

4 
3 
2 
1 



* 


negative 
i P+P 

* C+C 
■ Si+Si 

• Pb+Pb 

i 






— -f 


positive 
i P+P 

* C+C 
■ Si+Si 

• Pb+Pb 

" — "i— — 






• 

A^J-.;:- 


all charged 
1 P+P 

* C+C 
■ Si+Si 

• Pb+Pb 






FIG. 15: (color online) The scaled variance of the multi- 
plicity distribution of all charged particles produced in semi- 
peripheral Pb+Pb collisions as a function of the fraction of 
registered particles. The red diamonds correspond to the 
WA98 data [|. The error bars represent the systematic un- 
certainties, the statistical errors are smaller than the symbols. 
The straight line shows a linear fit to the points. See text for 
the details. 
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FIG. 14: (color online) The scaled variance of the multiplic- 
ity distribution of negatively (upper panel), positively (mid- 
dle panel) and all (bottom panel) charged particles produced 
in p+p, semi-central C+C, semi-central Si+Si and Pb+Pb 
collisions as a function of the fraction of nucleons which par- 
ticipate in the collision. The statistical errors are smaller than 
the symbols (except for the most peripheral points) . The hor- 
izontal bars indicate the systematic uncertainties (the upper 
one for the second point is off scale) . 



The NA49 detector registers a fraction p ~ 0.55 of all 
produced particles. The fraction p of registered particles 
can be decreased by cuts in the rapidity distribution. In 
the case of p+p interactions, the scaled variance of the 
multiplicity distribution was calculated in the following 
rapidity intervals Ay: (1.1, 2.6), (0.85, 3.8), (0.0, 3.8), 
(-0.4, 3.8) and (-2.4, 3.8) which correspond to fractions 
p of registered particles equal to 0.18, 0.25, 0.41, 0.48 
and 0.55, respectively. Instead of cutting in rapidity one 
can change the fraction of registered particles by select- 
ing a range in transverse momentum px of the registered 
particles. The scaled variance of the multiplicity distri- 
bution was then calculated in the following Apt (GeV) 
intervals: (0.0, 0.2), (0.0, 0.4), (0.0, 0.6) and in the full 
range corresponding to p equal 0.2, 0.37, 0.41 and 0.55, 
respectively. 



Fig. [TB] shows the scaled variance of the multiplicity 
distribution for negatively, positively and all charged par- 
ticles produced in p+p collisions at 158 A GeV as a func- 
tion of the fraction p of registered particles chosen by 
changing the rapidity and transverse momentum inter- 
vals. One can clearly see the linear dependence described 
by Eq. (p~5|) . The dependence of the scaled variance on 
the fraction of registered particles is independent of the 
choice of rapidity or transverse momentum as cut vari- 
able. 

The rise of the scaled variance of all charged particles 
with increasing acceptance can be presumably attributed 
to the effect of resonances which predominantly decay 
into oppositely charged particles. 

Multiplicity distributions were measured in p+p in- 
teractions in numerous experiments in a wide range of 
collision energies. An excellent compilation of p+p data 
can be found in [34], [3^, [3(| . Fig. [T7] shows the scaled 
variance of multiplicity distributions of all charged par- 
ticles in full phase-space as a function of the collision 
energy. The extrapolated result from NA49 (Eq. (fl5)) ) is 
also shown in Fig. [T7J The corrected scaled variance of 
all particles reaches a value of about 2 and fits well to 
the other data shown in this plot. 
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VII. COMPARISON WITH MODELS AND 
DISCUSSION 



Before closing our considerations we discuss possible 
origins of the strong dependence of Var(n) / (n) on the col- 
lision centrality. In this context, it should be stated that 
the strong centrality dependence of the scaled variance 
found in our experimental study was completely unex- 
pected. We do not know of any theoretical considerations 
that would have predicted the effect. However, after re- 
lease of our preliminary data [3j, l26j several attempts to 
explain them have been formulated 0, [H, H3, 0, [H| • 
Below the attempts are briefly presented. 

We start the discussion with a reminder that the multi- 
plicity fluctuations presented in this study were observed 
for a fixed number of the projectile participants. Al- 
though the average numbers of participants from the 
projectile and from the target are approximately equal 
in the collisions of identical nuclei, the number of target 
participants fluctuates even though the number of pro- 
jectile participants is kept fixed. We also note that the 
multiplicity fluctuations were measured in the forward 
hemisphere (1.1 < y c . m . < 2.6). Keeping these remarks 
in mind, we note that the large multiplicity fluctuations 
seen in the forward rapidity window are caused either by 
dynamical fluctuations of the particle production process, 
which are present even at fixed number of projectile par- 
ticipants, or the multiplicity fluctuations in the forward 
hemisphere result from the fluctuations of target partic- 
ipant number (the number is not experimentally fixed) 
which are transferred to the forward hemisphere due to 
an unknown collision mechanism. The mechanism has to 
strongly couple the forward and backward hemispheres. 

For the purpose of the further discussion, it is use- 
ful to divide the models of nucleus-nucleus collisions into 
three categories as proposed in [2(j. There are trans- 
parency models where the projectile (target) participants, 
which are excited in the course of interaction, mostly con- 
tribute to the particle production in the forward (back- 
ward) hemisphere. Furthermore, there are (rather un- 
realistic) reflection models when the projectile (target) 
participants mostly contribute to the backward (forward) 
hemisphere. Finally, there are mixing models in which the 
projectile and target participants contribute significantly 
to both the backward and forward hemispheres. 

In the transparency models, the observed multiplic- 
ity fluctuations can be caused only by the already men- 
tioned dynamical fluctuations of the production process. 
In the reflection and mixing models, the fluctuations of 
the target participant number contribute to the multi- 
plicity fluctuations in the forward hemisphere in addition 
to possible dynamical fluctuation. The models, which are 
represented in Figs. HU [H and [HI that is HIJING [27|, 
HSD [H|, and UrQMD all belong to the transparency 
class. The results of VENUS [3(| are very similar to those 
of HIJING, HSD and UrQMD, even though the correla- 
tion between the forward and backward hemispheres is 
somewhat stronger in VENUS. The model calculations 
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FIG. 16: (color online) The scaled variance of the multiplicity 
distribution of negatively, positively and all charged particles 
produced in p+p collisions as a function of the fraction of reg- 
istered particles chosen by change of rapidity (full symbols) 
or transverse momentum (open symbols) interval. The error 
bars represent the systematic uncertainties, the statistical er- 
rors are smaller than the symbols. The straight lines show 
the linear fits to the points. 



were performed with proper simulation of the VCAL re- 
sponse and the NA49 acceptance. The Ey etc , energy of a 
model event was calculated as the energy of the projec- 
tile spectators smeared by a Gaussian distribution with 
the width given by Eq. (jlOp . The models, which pro- 
duce approximately Poissonian multiplicity distributions 
independently of centrality, highly underestimate the ob- 
served multiplicity fluctuations in non-central collisions. 
The models are unable to reproduce even qualitatively 
the centrality dependence of the scaled variance. This 
remark applies not only to the Monte Carlo models as HI- 
JING 13, HSD [H, and UrQMD which are based 
on string excitation and decay, but to any transparency 
model which does not assume correlations among sec- 
ondary particles. For example, within a statistical model 
with fixed volume, where the electric charge is strictly 
conserved, the scaled variance of positive or negative par- 
ticles varies in the range 0.5 — 1.0 (3^|. The actual value 
depends on the thermal system's volume and the accep- 
tance in which the multiplicity fluctuations are observed. 

As already noted, the transparency models can pro- 
duce large multiplicity fluctuations, if there are built in 
correlations among produced particles. The problem can 
be formulated in a quite general way. 

The average multiplicity of secondaries is roughly pro- 
portional to the number of nucleons participating in the 
collision Np (cf . Fig. [5]) . Assuming that Np is propor- 
tional to the system's volume V, we have (N) — pV with 
p being the constant density of produced particles. Defin- 
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FIG. 17: The scaled variance of the multiplicity distribution 
of charged particles produced in p+p interactions as a func- 
tion of the center of mass energy. The full circles show a com- 
pilation of p+p results (cf. 34, 35, 36]). The open symbols 
show our results obtained by the extrapolation by changing 
the rapidity Ay or transverse momentum Apr intervals. 



ing the correlation function v(ri — r 2 ) = v(r) through the 
equation p 2 {ri,r 2 ) = p{ri) p{r 2 )(l + v{v\ - r 2 )), where 
P2(i"i,r 2 ) is the two-particle density, we get the desired 
formula 



Var(iV) 



= 1 + p / (fr v{y) , 



(16) 



which tells us that the multiplicity fluctuations are re- 
lated to the inter-particle correlations. The multiplicity 
distribution is Poissonian, if the particles are emitted in- 
dependently from each other (f(r) = 0). 

It is not difficult to invent a correlation function i>{r) 
which substituted in Eq. (|16|) reproduces the data shown 
in Fig. [IT] Some functions are discussed in [37] . The 
correlation function has to be positive at small distances 
(attractive interaction) and negative at larger distances 
(repulsive interaction). For Pb+Pb collisions the sign of 
the correlation changes at r sa 4 fm which corresponds to 
Np w 70 when Var(N)/(N) reaches its maximum. For 
r > (300) 1 / 3 fm « 7 fm the correlation function van- 
ishes. Possible mechanisms responsible for such a corre- 
lation function were discussed in [37]]. They include: a 
combination of attractive and repulsive interaction (the 
prediction of this model is shown in Fig. [18k). percola- 
tion, dipole-dipole interaction, and non-extensive ther- 
modynamics. 

Specific percolation models, which produce strong cor- 
relations among secondaries and thus provide large multi- 
plicity fluctuations, were discussed in [H, The color 
strings, which are stretched in the course of the interac- 
tion, form multi-string clusters. In p+p collisions there 
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FIG. 18: (color online) The scaled variance of the multiplic- 
ity distribution of negatively charged particles produced in 
Pb+Pb minimum bias collisions as a function of the number 
of projectile participants. The data are compared to the pre- 
dictions of the following models: (a) the model which assumes 
the inter-particle correlations caused by the combination of 
strong and electromagnetic interactions [321, (b) the percola- 
tion model [38I ]. (c) the transparency, mixing and reflection 
models [2Cj], see text for the details. 



are single-string clusters but in central Pb+Pb collisions 
the density of strings is so high that the strings over- 
lap and form one big cluster. The decrease of the scaled 
variance as a function of centrality is in this picture as- 
sociated with the percolation phase transition, i.e. the 
appearance of the large cluster. The prediction of the 
model [38| is shown in Fig. [TBb and reproduces the qual- 
itative behavior of our experimental results. 

While the transparency models require strong inter- 
particle correlations to comply with our data, the mix- 
ing and reflection models can produce large multiplicity 
fluctuations in the forward hemisphere due to the fluctu- 
ations of target participant number which is not fixed in 
our measurement. In Fig. 118b we show the prediction of 
the reflection and mixing models [13]. In the reflection 
model the target participants contribute to the forward 
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hemisphere, and thus the fluctuations of the target par- 
ticipant number produce the multiplicity fluctuations of 
secondaries. As seen in Fig. [18h . the reflection model 
strongly overestimates the observed scaled variance but 
the mixing models, where the target as well as projectile 
participants equally contribute to the forward and back- 
ward hemispheres, approximately agree with the data. 
The question arises whether the mixing models properly 
represent the dynamics of nucleus-nucleus collisions at 
high-energies. A recent theoretical analysis of d+Au col- 
lisions at RHIC |4(J suggests that the contribution of 
projectile participants indeed extends to the backward 
hemisphere while that of target participants to the for- 
ward hemisphere. 

In principle, the discussed models can predict multi- 
plicity fluctuations not only in Pb+Pb collisions but also 
in p+p, C+C and Si+Si reactions. It remains to be seen 
whether they can provide a consistent description of the 
fluctuations in all reactions. 

Finally we note that the transverse momentum fluctua- 
tions measured in nuclear collisions at 158 A GeV [1(3] and 
quantified by the show a similar centrality depen- 

dence with a maximum located close to that of the scaled 
variance of the multiplicity distribution. This behav- 
ior of the transverse momentum fluctuations as a func- 
tion of collision centrality was related in a superposition 
model [4l( to the centrality dependence of the multiplic- 
ity fluctuations using the correlation of average px and 
multiplicity observed in p+p collisions [Iol |. Moreover, 
we note that the percolation string model also describes 
qualitatively the centrality dependence of 3>(pt) [12] ■ 

VIII. SUMMARY AND CONCLUSIONS 

Multiplicity fluctuations in the forward hemisphere 
were analyzed by studying the scaled variance of the mul- 
tiplicity distributions in p+p, C+C, Si+Si, and Pb+Pb 
collisions. A strong centrality dependence is seen for neg- 



atively, positively and all charged particles. The number 
of projectile participants was used to determine the col- 
lision centrality. The scaled variance is close to unity for 
p+p, C+C, Si+Si and central Pb+Pb collisions. How- 
ever it increases significantly towards peripheral Pb+Pb 
collisions. The magnitude of the scaled variance is simi- 
lar for positively and negatively charged particles and is 
about 1.5 times larger for all charged particles. 

The string-hadronic models of nuclear collisions (with 
or without color exchange and without string fusion) 
predict no dependence of the scaled variance on the 
number of projectile participants, and thus the models 
qualitatively disagree with the data. A maximum of 
the scaled variance appears for semi-peripheral col- 
lisions when the forward and backward hemispheres 
are strongly coupled and fluctuations of the number of 
target participants contribute to the forward rapidity 
window. The observed scaled variance can also be 
reproduced within the models which assume strong 
dynamic correlations as, for example, in percolation 
phase transitions. 
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